The size of most massive neutron stars may reveal its exotic cores 
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The recent high precision observation of the most massive pulsar J1614-2230 with (1.97 ± 0.04) 
solar mass (Mq) was reported with a suggestion that many nuclear models which consider exotic 
particles in the core could be ruled out. However, many recent calculations could explain this star 
with various exotic particles, rendering the precise mass measurements insufficient to conclude on 
exotic cores. We examine the sensitivity of the radius of such a star to the details of its core. With 
our calculations and analysis, here we show that, for the most massive neutron star, with a precise 
observation of its radius it is possible to ascertain the presence of exotic cores. 



PACS numbers: 26.60.Kp, 26.60.-c, 97.60.Jd 
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The constituents of neutron star (NS), especially at its 
core remain as an intriguing subject for several decades. 
NS represents one of the final stage of the stellar evo- 
lution, and considered as a laboratory for the study of 
matter at super high densities. Observing the thermal 
spectra of NS results in a measurement of the gravita- 
tional redshift, which depends on both radius and mass 
of the star. With one observation the mass and radius 
can not be determined accurately However, precise 
mass measurements of NS have been possible Q , and sig- 
nificant work is in progress for measuring both mass and 
radius with more accuracy . Such measurements are 
crucial in improving our understanding of NS matter and 
the underlying interactions, because the observed masses 
and radii of NS constrain the equation of state (EoS) of 
NS or in general the nuclear matter and the finite nuclear 
properties such as neutron skin thickness in heavy nuclei 
Q, breathing mode resonance energies [1, Q etc. 

In an era where the nuclear models and their parame- 
ters are expected to account for the properties of infinite 
matter and neutron star d-lll], the new observation 
of the pulsar J1614-2230 with mass (1.97 ± 0.04) M© 
stands as a stringent test. This observation has a strong 
impact on the neutron star (NS) related studies mainly 
due to the ruling out of: (i) many of the standard and 
successful EoS by then (for e.g. PAL6 [H], FSUGold [|) 
which were yielding softer EoS and hence lower mass, 
and (ii) most of the EoS involving exotic matter such 
as kaon condensates or hyperons (for e.g. GM3 [l^, GSl 
[l6| ) which were "believed" to be in the core of NS. Soon, 
these claims were played down by J.M. Lattimer by 
stating "To rule out these exotic models fully you also 
need to know the radius of the star". He also added that 
it will be easy to tweak the parameters, to bring back 
the exotic particles and/or to reproduce the mass. This 
tweaking cannot be arbitrary due to several established 
constraints for the symmetry energy, EoS and finite nu- 
clear properties. Many of the recent work were in this 
direction where few more EoS were introduced jl8| and 
in several models the parameters were adjusted [3, fioj . 
to reproduce the 2Mq NS. However, in the context of 
confirming/ruling-out exotic cores in neutron stars, the 
precision required in observation of radius of most mas- 



sive neutron stars, is not reported so far. 

In our recent work [l^ we found that for a NS with 
maximum mass 2Mq, without antikaons (G2) we get the 
radius is 11.03 km whereas, with K~ (FSU2.1) and both 
antikaons (Gl) the radii are 11.42 km and 12.55 km, 
respectively. This sensitivity of the radius to the pres- 
ence of exotic cores motivated us for a theoretical survey 
which may demand more precise observation of radius 
along with the mass of such a heavy NS. In the follow- 
ing text we outline the results from (i) our calculations 
with recent versions of extended rclativistic mean field 
(RMF) models which are one among the most reliable 
models for NS EoS, (ii) variety of interactions in liter- 
ature and (iii) a fiducial model with parabolic EoS, all 
successfully explaining a static 2Mq NS. In this work 
we have assumed that the maximum in observed neutron 
star masses [(1.97±0.04)Mq] to be same as the maximum 
mass allowed by a chosen EoS. It is likely that latter can 
be at least a few tenths of solar mass larger Impli- 
cations of this on our results are studied with the fiducial 
model. 

In the RMF model, we have chosen few representa- 
tions, which signify as extrema with and without an- 
tikaons {K^,K^), that can yield 2Mq as the maximum 
mass of NS. The model Lagrangian and details of calcu- 
lation are explicitly explained elsewhere [ll|, [H, [2l| . Us- 
ing the Lagrangian density, we obtain the energy density 
and pressure (EoS) with or without antikaons. Once the 
EoS is defined, we use the Tolman-Oppenheimcr-Volkoff 
equations [l^ll^l to get the mass-radius relation of a NS. 

The mass-radius relation, thus obtained, using the pa- 
rameters NL3[l3l, Gl, G2 [nj, and FSU2.l[i3 are shown 
in Fig. 1. Among these chosen parameters, G2 yields 
the softest EoS. In general, inclusion of exotic particles 
soften the EoS and a softer EoS results in lesser mass for 
NS. Without inclusion of any exotic particle, G2 yields 
a mass consistent with the 2Mq NS. Thus G2 represents 
one extremum where exotic particles cannot be accom- 
modated. The other extremum should correspond to an 
EoS with reasonably extreme value of Uk (if we consider 
only antikaons as the exotic particles) yielding a 2Mq NS. 
NL3 represents this extremum and its EoS is the stiffest 
among the chosen parameters. The other parameters Gl 
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FIG. 1: The mass-radius relation from RMF models [l^. Different curves represent the calculations using different kaon optical 
potentials (Uk) and with different parameter sets. For each parameter set, solid black line represents nucleonic (non-kaonic) 
phase, lines with different patterns and colors represent the phase with K~ and the corresponding small dotted lines represent 
the phase with both antikaons {K~ ,K^). The different patterns and colors represent the strength of the kaon optical potential 
Uk {\Uk\ quantifies the influence of kaons) as specified in the inset . The sensitivity of mass-radius relation to the parameter U k 
depends on the stiffness of symmetry energy. For e.g., FSU2.1 has softer symmetry energy and hence it is not much sensitive 
to Uk- The onset of kaon condensation is governed by the interplay between density dependance of, symmetry energy and 
EoS [U]. The delayed onset of K~ in FSU2.1, due to the above reason, forbids the onset of R'^ . The solid circles represent 
the maximum mass in every case. Mass is given in units of solar mass (Mq). The shaded region correspond to the recent 
observation of (1.97 ± O.O4)M0 neutron star. 



TABLE I: Radii of neutron star of maximum mass (1.97 ± O.O4)M0 resulting from different models/parameters and with 
different compositions as reported in the recent literature. 



Model/Parameter 


Radius (km) 


Composition 


APR (NS) [22] 


11.0 




G2 [21] 


11.03 


n,p, e",/i" 


lU-FSU [9] 


11.2 


n,p,e",/i" 


FSU2.1 [19| 


11.42 


n,p,e",^", A"" 


Gl [21J 


12.45 


n,p,e",/i",/s:" 


Gl [19| 


12.55 


n,p,e~ ,K- ,K0 


NL3 (This work) 


13.04 


n,p,e~,fj,~,K-,K° 


GlendnhS-fFth [23| 


11.48 


n, p, e~ , /i~ +Hyperons 


DD2H[24| 


11.6 


n,p,e~ , -f Hyperons 


QMC700/NY[25l 


12.5 


n,p,e^ , fi^ +Hyperons 


QMC-HF NYfi:7-7r-m; [26| 


12.38 


n,p,e~ , /i^ , 7r-|-Hyperons 


QCD (a4=.62 and 77 = -1)[27| 


10.60 


n,p,e~ , fi~ , K~ + Quarks 


APR (HS) [22| 


12.2 


n,p,e~,fi~ + Quarks 


DD2H-NJL, set A[24] 


11.8 


n,p,e~ , +Hyperons-|- Quarks 


DD2H-NJL, set B[24| 


12.15 


n,p, e~, /i^-fHyperons-|- Quarks 


GMl-fPNJL[28] 


13.2 


n, p, e~, /i~+Hyperons 4- Quarks 



and FSU2.1 fall in between the extrema and have similar 
EoS but significantly different symmetry energy contri- 
bution. From Fig. 1, one can observe a general trend that 
the radius (i?) corresponding to maximum mass is lower 
for a softer EoS. However, this scenario can change in 
the presence of antikaons. Though the antikaons soften 
the EoS, they can increase or decrease the central baryon 
density (pc) [2l| and hence R {R (x 1/pc)- This results 
in, different patterns of variation of i?, as shown in Fig. 1. 

These results suggest that a 2Mq NS can be explained 



with different quantity of antikaons, and the radius {R2) 
corresponding to a NS with a maximum mass of 2Mq is 
different with different parameter sets. Our results from 
RMF models show that if R2 > 11.04 km then the exotic 
core must be present. In our calculations exotic core 
is restricted to antikaons, but these calculations should 
be extended to include hyperons and quarks. Similar 
calculations for the 2Mq NS with other models, and a 
careful scrutiny of the effect of rotation on massive NS, 
are needed to have a complete outlook. 
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FIG. 2: The mass-radius relation from fiducial model, (a): Schematic representation of EoS from our fiducial model comprising 
one or two parabolas represented by dashed and solid lines respectively. The dashed line does not change abruptly as a function 
of density. This could effectively mimic an EoS of NS with nucleonic phase only, with a single parameter (a) representing 
the stiffness of the EoS. In case of EoS with two parabolas with slopes a\ and 02, the point of intersection (eo,Po) represents 
the onset of exotic phase whence the EoS changes abruptly. Both these EoS can yield same mass but different radii for a 
neutron star, (b): The radius of the neutron star (when its mass is maximum) plotted as a function of the slope of parabola 
when the EoS is assumed to be a single parabola. Green region is where the EoS yield a neutron star with maximum mass 
(1.97±0.04) A'/q , and blue squares indicate few selected RMF models as shown in the labels. The first and second stars represent 
the cases where the maximum masses are 2.2Mq and 2.4M0, respectively, (c): Sample results for the radius of a neutron star 
with maximum mass (1.97 ± O.O4)M0 obtained from an EoS with two parabolas, plotted as a function of the slope of second 
parabola (02) representing the exotic phase and the point at which the second parabola starts to contribute (eo)- 



In Table HI we have compiled R2 calculated with differ- 
ent interactions and compositions as quoted in recent lit- 
erature. We have excluded the cases which yield the max- 
imum mass higher than 2.01Mq. For further discussions, 
it is convenient to define two radii R2 and (in km) 
for the cases of NS without exotic core and with exotic 
core respectively, corresponding to a maximum mass of 
(1.97±0.04)Mq. From TableUwe have 10.6 < i?f < 13.2 
and 11.0 < R2 < 11.2. Thus we can claim that if R2 
were observed beyond the range of R2 then the exotic 
core must be present. Otherwise, the existence of exotic 
cores will remain ambiguous. To strengthen this claim, 
one has to look more into the systcmatics, beyond relying 
on few theoretical results available till date. 

As a first step, we start with a fiducial model where 
the EoS is a single parabola represented by the dashed 
line in Fig 2a. The resulting radii are shown in Fig. 2b as 
a function of slope of the parabola (a) . Smaller the value 
a, softer the EoS, more is the central density and hence 
we obtain a smaller NS. By varying a, we could fit EoS 
of different RMF models. The radii obtained with such 
fitted EoS are quite consistent with those obtained from 
the actual EoS, justifying our approach. These radii are 
marked in Fig 2b. For a NS with mass (1.97 ± 0.04)AfQ, 
our fiducial model calculations predict that 2.7 < a < 
2.93 and correspondingly 10.87 < R2 < 11.27. However, 
this range is very inclusive because many values of the 
parameter a yield quite unrealistic EoS which will not 
satisfy several experimental constraints. 

In our next step we consider two parabolas for the EoS 



TABLE II; Radius of a neutron star of mass (1.97±O.O4)M0, 
obtained from different models assuming this mass as the 
maximum mass. 

Model Radius (km) 

Without exotic core With exotic core 
RMF (Fig. 1) ILOS 11.43 <R< 13.04 

Parabolic EoS 10.87 <R< 11.27 9.79 <R< 13.22 
FromTable[T] 11.0 < < 11.2 10.6 < i? < 13.2 



represented by the solid line in Fig 2a. Sample results for 
R2 with a fixed slope of the first parabola (oi) are given 
as a function of slope of the second parabola (02), and the 
energy density (sq) at which the second parabola starts 
to contribute to the EoS. We have considered a broader 
range of Eq whereas the lower limit of Eq shall correspond 
to quarks which typically start to contribute from ~ 100 
MeVfm"^. £0 and 02 varied in tandem can effectively 
represent the parameters of exotic matter EoS. In case 
of antikaons we found that they simulate the variation in 
Uk- The important feature of Fig. 2c is that the radius 
is quite sensitive to the parameters of the exotic matter 
EoS and the variation in Eq and 02 leads to a change in 
radius up to ~ 2.5 km. We can see that for a given ai 
and 02, Rf increases with eq. Lower values of Rf are 
due to the early onset of exotic phase. The range of R2 
obtained for the chosen ai value is presented in Fig. 2c 
but a complete range for all allowed values of oi, 02 and 
£0 is given in Table |lll where we summarize all of our 
results. 
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Our fiducial model suggests that without exotic core 
we have a narrow allowed range of i?2- If R2 were ob- 
served beyond this range, for e.g. R2 ~ 12.3 ± 1 km or 
10.3 ± 0.5 km, then an exotic core must be present. The 
upper limit of this allowed range could be higher, if (i) a 
NS is observed with mass above 1.97 ± .G4M0 or (ii) we 
consider that the maximum mass given by the EoS has to 
be higher than the maximum in observed mass. Any in- 
crease in mass will result in the increase of radius without 
the exotic cores [See Fig. 2b] . Consequently, cases (i) and 
(ii) will have no bearing in the lower limit of the range, 
i.e., the least radius that could be supported by a NS 
without exotic cores. Hence, it is clear that if i? < 10.87 
km then the NS must have an exotic core. The range of 
radii obtained by the fiducial model may vary, if one con- 
siders two polytropes instead of two parabolas but with 
two more parameters to adjust. However, from Table HIl 
we clearly see that the ranges of R2 given by our fiducial 
model comprise the ranges obtained from our RMF calcu- 
lations and other recent results. Hence safer conclusions 
could be drawn from these broader ranges of i?2- If we 
consider more realistic form for the EoS and incorporate 
the experimental constrains, for the saturation proper- 
ties defining the EoS and the high density behaviour of 
the EoS (For e.g. [3lj), the range of i?^ may become 



narrow. 

For a static neutron star with maximum mass (1.97 ± 
O.O4)M0, we conclude that a pure nucleonic star can 
have only a narrow range of radius (For e.g. 10.8 km 
< i? < 11.3 km) and the presence of exotic cores widen 
this range. Beyond the above-mentioned range, exotic 
cores must be present. If the maximum mass given by 
an EoS has to be larger than the maximum in observed 
NS masses then the upper limit of above range could be 
different but the constraint on the lower limit remains 
intact. More rigorous calculations are needed to obtain 
precisely the range of radius allowed for the most massive 
neutron star without exotic particles. 

The 2Mq NS does not rule out exotic cores and hence 
the constituents of NS still remain a puzzle. To over- 
come this, more observations and calculations should fo- 
cus on the radius of such massive neutron stars. A pre- 
cise measurement of mass and radius of a lower mass 
neutron star will not be useful confirm the presence of 
exotic cores. Our calculations and analysis predict that 
a precise observation of radius of a 2Mq NS is capable of 
confirming the presence of exotic cores. This in turn can 
shed light on the high density behaviour of matter and 
the phase transitions that can happen at such extreme 
environments. 
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